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Abstract 

We study the isospin dynamics in fragment formation within the framework of an analytical 
model based on the spinodal decomposition scenario. We calculate the probability to obtain frag¬ 
ments with given charge and neutron number, focussing on the derivation of the width of the 
isotopic distributions. Within our approach this is determined by the dispersion of N/Z among 
the leading unstable modes, due to the competition between Coulomb and symmetry energy effects, 
and by isovector-like fluctuations present in the matter that undergoes the spinodal decomposi¬ 
tion. Hence the widths exhibit a clear dependence on the properties of the Equation of State. By 
comparing two systems with different values of the charge asymmetry we find that the isotopic 
distributions reproduce an isoscaling relationship. 
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I. INTRODUCTION 


In the last years a widespread attention has been devoted to the role played by the 
isospin degree of freedom in the heavy-ion reaction physics. The interest on this subject is 
twofold: the knowledge of the symmetry term in the Equation of State (EOS) of asymmetric 
nuclear matter, which is a fundamental ingredient in astrophysical investigations [lj, and 
the thermostatistical properties both at equilibrium and out of equilibrium of systems with 


two strongly interacting components 


Both the interests concern systems 


faraway from the physical conditions of ordinary nuclear matter. 

Thanks to the availability of high-performance 47r-detectors for the investigations of 
heavy-ion collisions at intermediate energy BEQ1I3, recent experimental results can 
provide new insights about isospin effects on the nuclear dynamics. In particular, for mul¬ 
tifragmentation processes we can obtain information about highly excited two-component 
systems and their subsequent decomposition. Statistical models have been extensively ap¬ 
plied to the description of experimental data, also for isospin observables [13], and some 
conclusions have been drawn on the behavior of charge asymmetric systems. These mod¬ 
els, however, imply the achievement of the statistical equilibrium for the nuclear system. 
Then, it would be highly desireable to have some insight on the path followed by the system 
to attain equilibrium, if this occurs. Further, it would be of great advantage to envisage 
some observable, which preserves memory of the dynamical processes occurred during the 
fragmentation. 

In this paper we present an analytical description of the disassembly of excited nuclear 
systems formed during the collision of heavy ions, in terms of the occurrence of nuclear 
matter instabilities. Our approach accounts for the source of the density fluctuations oc¬ 
curring when the system enters the spinodal instability region of the density-temperature 
phase diagram, and describes the growth of the fluctuations with time until they cause the 
decomposition of the system. This approach is a generalization to include the isospin degree 


of freedom, of the model developed in Refs. [1J, |15|] for symmetric nuclear matter basically. 


This gives rise to a substantial improvement of the model, with new valuable results. Such 
extension allows us to investigate separately fluctuations of the neutron and proton densi¬ 


ties and their interplay. Following the procedure introduced in Ref. |14j, we identify the 


pattern of the domains containing correlated density fluctuations, with the fragmentation 
pattern, and can make predictions on the isotopic distributions of the fragments. Moreover, 
we include in the present treatment the Coulomb force according to the approach outlined 
in Ref. [16]. Its effects on the isotopic distributions turn out to be sizeable. 

Our results essentially refer to the distributions of the fragments just after the early break¬ 
up of the system. So our approach can be considered complementary to dynamical model 
calculations based upon semiclassical kinetic equations for one-body phase-space density, 
(for a review on dynamical models see, e.g., Refs. 17, 0, 519]), as far as the description 
of the early fragmentation mechanism is concerned. The advantage here is that one can 
make signihcant predictions on observables of experimental interest on an analytical basis. 
This allows us to directly relate the results obtained to the EOS properties and the features 
of the spinodal mechanism. In our scheme the onset and the growth of the fluctuations 
about the mean phase-space density in unstable situations, are self-consinstently treated. 
The self-consistency condition is provided by the fluctuation-dissipation theorem. Whereas 
all the processes, which take place before the system enters the spinodal instability region 
and after the break-up, are beyond our approach. Therefore the mean values of density, 
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temperature and asymmetry of the nuclear medium when the system starts to break up are 
taken from calculations performed within dynamical models. 

On the other side, a dynamical model, which appropriately incorporates the effects of the 
fluctuations, might give a detailed description of the whole history of a collision between 
heavy ions. Therefore, it can be of interest to compare the results of our approach with 
those obtained by numerical solutions of microscopic transport equations, also to connect the 
results of the simulations to what is expected in a pure spinodal decomposition scenario. The 
comparison will be done with the isotopic distributions for the primary fragments, calculated 
in the dynamical stochastic mean-held (SMF) approach of Refs. [hj,|20(. In particular, we will 
consider the ratio, for a given value of the proton number, between the isotope yields from 
two different reactions. This quantity represents a straightforward mean to compare isotopic 
distributions, since it is experimentally found to obey a simple relationship (isoscaling), as 
a function of the proton number and neutron number |9, llOJ, |2jJ, u 
the dependence of the isoscaling parameters on the EOS considered. 

In Sec. II we outline the extension of the formalism developed in Ref. 
isoscalar density fluctuations, to include the isospin degree of freedom. In Sec. 


We will also discuss 


14] only for 

III we 

discuss the results of our calculations and their comparison with the calculations performed 
in Ref. 23|| within the SMF approach. Finally, in Sec. IV a brief summary and conclusions 
are given. 


II. FORMALISM 

A. Time evolution of density fluctuations 

We study the density fluctuations by introducing a self-consistent stochastic field acting 
on the constituents of the system. The time evolution of the fluctuations is described 
by a kinetic equation, within a linear approximation for the stochastic field. The growth 
of fluctuations is essentially dominated by the unstable mean field. Thus we focus our 
attention on the behavior of the mean field and neglect the collision term in the kinetic 
equation. Collisions would mainly add a damping to the growth rate of the fluctuations and 
should not change the main results of our calculations, at least at a qualitative level. 

The additional stochastic mean field, which we assume having a vanishing mean, will in¬ 
duce fluctuations of the proton and neutron densities, £pi(r, t ), with respect to their uniform 
mean values Qi { i — 1,2 for protons and neutrons respectively ). We assume that at the 
time t = 0, given density fluctuations 5gj{r,t = 0) are present in the system. The equations 
for the Fourier coefficients of <5pj(r, t) for t > 0 a re g iven by a generalization of the equation 
for the isoscalar density fluctuations of Ref. [lj, |24j. They read 

<5pi(k, t ) = <5pj(k, t = 0) — Ik i <5p/(k, t = 0 )D^(k, to = 0) / A j(k, t — t ') dt' 

J ’ Jo 

+A f'D id {k,t,^dWjik, If). (1) 

Jo 

where the 2x2 matrix in the isospin space, D lJ (k, t — t '), is the density-density response 
function and Dij(k,to) its time Fourier transform. For symmetry reasons the response func¬ 
tion and its Fourier transform depend only on the magnitude of the wave vector. In the 
last integral dWj(k,t') gives the contribution of the j-component of the stochastic field in 
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the interval dt' . Since the stochastic field is real IF*(k, t) = Wj(—k, t). The real and imagi¬ 
nary parts of the Fourier coefficients Wj(k, t) are indipendent components of a multivariate 
stochastic process [25|, with 

< [* dWi(k, If) f dWA-kfl") >= f dlfdf BijfclfX) (2) 

Jo Jo Jo 

defining the correlator for the stochastic field. Angular brackets denote ensemble averaging. 

In the mean-field approximation the response function obeys the following set of equa¬ 
tions 

D id {k,u) = Di°\k,u)Sij + S iD^ik^AiMDui^u) (3) 

where D[°\k, u) is the non-interacting particle-hole propagator and A% } i{k ) is the Fourier 
transform of the nucleon-nucleon effective interaction. 

In Ref. 0 it has been shown that, in the case of isoscalar fluctuations in symmetric 
nuclear matter, a white-noise hypothesis for the stochastic field can be retained for values 
of temperature and density sufficiently close to the borders of the spinodal region. In such 
situations the imaginary part of the response function displays a sharp peak dominating the 
particle-hole background at a value of u <C kvF- This is due to the occurrence of a pole on 
the imaginary axis of u, that corresponds to isoscalar fluctuations, at a distance from the 
origin that is much smaller than the values of kvp- The position of this pole determines the 
time scale characteristic of the response function. 

However, when one wants to investigate the properties of neutron and proton distribu¬ 
tions, as we do in the present study, one should consider also the effects due to the isovector 
fluctuations. Even though isoscalar modes are the dominant ones, since they are unstable, 
isovector fluctuations contribute to the width of the isotopic distributions of the fragments 
formed in the spinodal decomposition process. In asymmetric nuclear matter isovector and 
isoscalar fluctuations are coupled. However one can still distinguish oscillations with neu¬ 
trons and protons moving in phase (isoscalar-like) or out of phase (isovector-like). Let us 
first concentrate on the properties of the isoscalar-like modes. 


1. Isoscalar-like fluctuations 

The position of the pole uj = for the unstable isoscalar-like mode is given by the 
imaginary root of the equation 


det|(5jj — Df\k, u)Ai t j{k) \ = 0 


(4) 


The quantity T/ c is the damping or growth rate (depending on its sign) of the density 
fluctuations. In evaluating it, we use the expression of Df\k, u) for u <C kvF HI 


Dr{k,uj) ~ 


d,a, - F(/J % ’ 


where the effective chemical potential fi, of neutrons or protons is measured with respect to 
the uniform mean field U t (o \, Q 2 ) of the unperturbed initial state and F(f3fli) is the function 


F{m = 


e PF + 1 ’ 
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with f3 = 1/T being the inverse temperature (we use units such that h — c = ks — 1). 

Substituting into Eq. © the response function Dij(k,t — t') calculated with these ap¬ 
proximations, the equation for the fluctuations 8gi(k,t) becomes 

Sgi(k,t) = 5pi(k,t = 0) + T, jt iC i) i(k)Df 1 (k,u = 0)<%(k,t = 0)-^-(e rfe * - 1) 

’ ' 1 fc 

+Ej Cij(k)e rkt f e~ Tkt ' dWj(k, t '), (5) 

J 0 

where Cij(k) are the residues, times (— i), of the components of the response function at 
the pole u = iT*,. They have the relevant property 


det|Cij(A;)| = 0. 


( 6 ) 


The explicit expression of the inverse of the response function for ui = 0 is 


D a }(k,u = 0) = — + Aij(k) 


h3 


djlj 


dpi 


For isoscalar-like fluctuations Wj(k,t') represents a Gaussian white noise 0 . The prob¬ 
ability distribution of density fluctuations, P[8gi(k, £)], is given by a product of Gaussian 
distributions. Each single factor corresponds to the stochastic process of Eq. © for a given 
wave number k jlil fidl]. with the covariance matrix 


cr 


lj(k,t) = T, l)m Ci }l (k)B l)in (k, t)C m j(k)-^—(e 2Vkt - l) . 

Z1 k 


(7) 


For simplicity, we have assumed that the initial fluctuations are negligible afj(k,t) ~ 0. 
Whenever it is necessary, a nonvanishing covariance can be easily introduced. 

The probability distribution P[8gi(k, £)] is completely determined once the covariance 
matrix crfj(k,t) is known. According to the procedure usually followed when treating in¬ 
stabilities by exploiting the fluctuation-dissipation theorem, see e.g. Refs. Q, 27|, we 
determine the coefficients B it j(k,t) as functions of 0 1; g 2 and T for the system at equilib¬ 
rium, then we extend the expressions so found to non-equilibrium cases. Since the relevant 
values of the wave vector k turn out to be such that the quantity kvp is of the same order 
of magnitude as T, the limit uj/kvp <C 1 also implies co/T 1. In such case, the clas¬ 
sical limit u/T -C 1 (or \T),{t)\/T -C 1) can be taken when evaluating both sides of the 
fluctuation-dissipation relation. Then, we get 


< <5ft(k,i)5^(-k,i') >= -TD itj (k,t~ t'). (8) 

The equation for the equilibrium fluctuations can be obtained from Eq. © by shifting the 
initial time t = 0 to — oo. By exploiting Eq. (JHJ) we can obtain the following relation between 
the coefficients Bij(k,t) and the functions Cij(k): 


Ei,mCi,[(k)Bi ym (k,t)C m j(k) = -2 TCij(k). 


(9) 


From this equation we can see that B it j are constant and depend only on the magnitude k 
of the wave vector, as it is expected for symmetry reasons. Following Refs. [261,127( (see 
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also the discussion in Ref. 3| on this point) we assume that the relation © is valid also 
in instability situations. In such a way, the covariance matrix © acquires the form 


cr 


■i(M) = - l) , 


( 10 ) 


and is completely determined both for stable and unstable situations. We notice that, for 
the isoscalar-like mode, crf 2 (A;) = a 2 i(k) i s positive. In fact proton and neutron densities 
oscillate in phase, although with different amplitudes in general. However, the ratio between 
amplitudes, of \(k)/erf 2 (k), is found to be larger than the initial proton to neutron ratio, 
thus leading to the formation of more symmetric fragments, the so-called isospin distillation 
effect jd]. 


2. Isovector-like fluctuations 


Now we turn to consider the isovector-like modes. In this case the frequency of the modes, 
uj™ is real, i.e. we have stationary oscillations. The position of the pole is given by the other 
solution of Eq. ©• However, we add a small negative imaginary part — T™ to the position 
of the pole, taking into account that here we are neglecting nucleon-nucleon collisions and 
finite size effects. Correspondingly the imaginary part of the response function acquires the 
width Yfl. 

The contribution of isovector-like fluctuations to the covariance matrix can be 

written as it follows: 


cr. 


hJ 


(M) =4 ^ m Ci}(k)CZ/k)e~ 2r ^ 

x f dtidt[ \e r ™ (t 1 +t>1 ^B™ m (k, ti,t[) sin (2 u l k v (t - H)) sin (2 J%{t - tf))] , (11) 


where C\^{k) are the residues at the pole and denote the contributions from 

the isovector-like fluctuations to the stochastic held. 

To determine the amplitude of the stochastic held we essentially follow again the deriva¬ 
tion presented above. By exploiting the fluctuation-dissipation theorem, now in the limit 
u/T » 1 (since the frequency of the isovector vibrations is rather large with respect to the 
relevant values of T), we obtain for values of to close to the pole the relation: 




IV 

l.m 


(k ,<-H) = 2 rfq" 


2(U 


(w-<4’') 2 +(I?) 


( 12 ) 


where we have added a Lorentzian factor to the right hand side in order to restrict to a small 
region about c v™ the contribution from the isovector-likc pole to the time Fourier transform 
of Bi v m (k, to). In this way the correlator B\ v m { k, t\ — t \) for the stochastic held results to be 
proportional to e -I 4b tl ~ i il. This means that the isovector-like stochastic held is given by a 
coloured noise, at variance with the isoscalar case. 

Substituting the time Fourier transform of Eq. m into Eq. (ED, and retaining only the 
leading term of the expansion in powers of (we obtain for the covariance matrix 
the expression 


<T?j(k,t)=C% 


1 — e 


-2ir t 


-2rr te 


-2 n v t 


+ 0(ft) 2 ), 


/ ' 
uk 


(13) 
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whose asymptotic value is given by 

^■(*) = ££■(*)■ (14) 

We notice that, for isovector-like fluctuations, crf 2 (k) = cr| 1 (/c) is negative. Indeed neutron 
and proton densities oscillate out of phase. 

The covariance matrix of Eq. m refers to equilibrium fluctuations at given values of 
density and charge asymmetry. It can be directly obtained by means of the fluctuation- 
dissipation relation in the case of a purely real pole ( T™ —s► 0 ). 

We finally remark that the covariance matrix of Eq. m is obtained in the limit T —► 0 
and, in addition, it does not depend on the width T™ of the isovector-likc resonance. This 
implies that the density fluctuations of isovector-like nature, we are considering, have a 
quantum origin. 


B. Size distributions 

Now we describe the procedure to determine the distribution for the size of the correlation 
domains. We closely follow the derivation given in Ref. □ for isoscalar density fluctuations, 
and we limit ourselves to outline the steps relevant to the present more general treatment. 
We distinguish the fluctuations of the proton density from those of the neutron density. 


1. Correlation lengths 


The probability distribution for the sizes of the domains where the fluctuations are cor¬ 
related, b\ and b 2 for protons and neutrons respectively, can be obtained by means of the 
functional integral 

P(b 1 ,b 2 ,t) = I d[5g i (r,t)]6(b 1 - J drdr'dg^r^^fi^Sg^r^^fiir')^ 

d(b 2 - J drdr , Sg2(r,t)f2< y r)Sg2(r , ,t)f 2 (r')^P[6gi(r, t)} , (15) 

where P[5gi(r,t)] is the probability distribution for the density fluctuations and /)(r) are 
suitable weight functions. Moreover, we assume that the dynamical correlation lengths for 
proton and neutron density fluctuations, < bi > and < b 2 >, coincide 

m = / ^d,i(M)i/i(*or= / (io 


where fi(k) are the Fourier transforms of the weight functions. In this way we assume that, 
on average, neutrons and protons are correlated within the same domain. We will see in the 
following how this can be related to the average isospin distillation effect in the formation 
of fragments. 

Following the procedure used in Ref. 

P(bi,b 2 ,t ) the equation 

P(bi, b 2 , t) — 77-T7IV n- 1 - T r^~~ ' ex P 1 ~ 


14l | we obtain for the probability distribution 
[bi + b 2 \' 


2vr L{t) [&i + b 2 \ 

1 fa - b 2 f 

4L(t)7(t) [bi + b 2 ] 


4L(f) 


xexp — 


(17) 


7 










(18) 


where the parameter 7 (t) is given by 


/ dkaj^jk, t)\fi(k)\ 2 J dka’ 2 2 (k, t)\f 2 (k)\ 2 
f dkaf 1 (k,t)\f 1 (k)j 2 f dkal 2 (k,t)\f 2 (k)\ 2 ' 


At variance with the case of isoscalar fluctuations, the distribution P(b 1 , b 2 , t) depends on the 
weight functions fi(k). These functions, to some extent, are arbitrary, the only requirement 
is that the integrals containing them should converge. For simplicity, we assume \fi(k)\ 2 = 
di\f(k)\ 2 . For the functional form of \f{k)\ 2 we choose the simplest one: \f(k)\ 2 = 1/k 2 . 
This choice is also supported by the fact that for equilibrium fluctuations the integral of 
the variance weighted with 1/k 2 gives the correct value of the correlation length jlj]. In 
addition, we have found that for the physical situations considered in this paper, the value 
of the parameter 7 (t) to a large extent is insensible to the particular form of the weight 
function \f(k)\ 2 . 

From the probability distribution of the domain sizes we can obtain the distribution of the 
numbers of correlated protons Z and neutrons N , assuming the correlation domains to be 
spherical. The relations between Z and b\, and N and b 2 can be expressed as 61 = 2r m Z 1 / 3 
and b 2 = 2r 02 N 1//3 , where r 0i is the mean interparticle spacing for nucleons of the z-species, 
calculated at the actual values of asymmetry and density (when fragments are formed), that 
are different from asymmetry and density of the initial matter. The fact that the fragment 
size is related to the correlation length can be considered as a reasonable assumption in 
situations where isoscalar-like modes are the dominant ones, as in fragmentation processes. 

So, since on average &i is equal to b 2 , we obtain: r 01 /r 02 = (p 2 /pi ) 1//3 =< N 1//3 > 
/ < z 1 / 3 >, where pi are the densities calculated at the time fragments are formed. In 
this way the ratio r 0 i/r 02 can be related to the average asymmetry of the liquid (fragment) 
phase, obtained after the distillation process has occurred. One can consider, for instance, as 
average fragment asymmetry, values extracted from dynamical SMF simulations for primary 
fragments ) 8 ]. 

Then, the probability distribution of Z protons and N neutrons contained in a correlation 
domain, acquires the form 


P(Z, N,t) = 


1 rg AiA 2 1 1 

9 vUfy [AxZ 1 / 3 + A 2 iW 3 j JzNfl - 3 y 7 (t) eXP 

f r 0 1 [A 1 Z 1 / 3 — X 2 N^ 3 ] 2 \ 

xexp {-mm Puzos + a 2 at/3] j 


_n)_ 

2 L{t) 


+ a 2 a^ 1 ^ 3 ]) 

(19) 


with A i = roi/ro, where ro is the mean interparticle spacing for nucleons of both species. 


2. Correlation volumes 


One may also assume that the size of fragments is directly related to a correlation volume 
V, instead of a correlation length. Equation (El) can be rewritten for the correlation volumes, 
just replacing bi and b 2 with V\ and V 2 . Then the probability distribution, after some algebra, 
reads: 


P(Z, N, t) = 


2tt V(t) \p 2 Z + Pl N] y7^) eXP ( 4 V(t ) [Z/Pl + " /P2] ) 

xex 7_Z_ 1 1gMjyyfajy 

P V 4V ( t) 7 (t) [Z/pi + N/p 2 \ J 


( 20 ) 
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where V is the average correlation volume for nucleons of both species. For not too large 
asymmetries, this can be rewritten in the following form: 


P(Z, N,t) = 


ttAA 
xexpl — 


7(f) 
A 


=exp — 


A 
2 A 


1 r N 


2 A 7 (t)L A 


Z i2 

- a 


( 21 ) 


where a = (p 2 — Pi)/(P 2 + Pi) represents the average asymmetry of fragments and A is the 
average mass. 


III. RESULTS 

In our calculations we have adopted a schematic Skyrme-like effective interaction, that 
can be expressed as a sum of two terms 


Ai,j(k ) — A(k) + Sij(k). 


For the symmetric term A[k) we use the finite-range effective interaction introduced in Ref. 

0 : 

^(*) = (/i2- + ( ff + i)2L^) e - 3 ^ i (22 ) 

V Qeq Peg 


with g = Qi + g 2 and 


A = -356.8 MeV, B = 303.9 MeV, a = 


These values reproduce the binding energy (15.75 MeV) of symmetric nuclear matter at 
saturation (g eq = 0.16 fm -3 ) and give an incompressibility modulus of 201 MeV. The width 
of the Gaussian in Eq. (1221 1 has been chosen in order to reproduce the surface-energy term 
as prescribed in Ref. |29|. 

The isospin-dependent part, S tJ (k), contains three different terms 


s ^- d ^ + ^ Dk2 ^v cm ,,, 


(23) 


with T\ — 1 and r 2 = —1. The double derivative of the potential part of the symmetry 
energy density, £ sy mm , is calculated in the unperturbed initial state. For the coefficient of 
the isovector surface term we use the value D = 40 MeV • fm 5 • Concerning the Coulomb 
interaction, a mean-held exchange contribution 


= -aW 


/ 3 \ V 3 2 - 2/3 
e 0 ! 


is added to the bare Coulomb force. 

In order to stress the effects of the asymmetry of the nuclear medium, we will present 
results obtained with two different parametrizations of the symmetry energy: one with a 
stronger density dependence ( “superstiff” asymmetry term ) and the other one with a weaker 
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FIG. 1: The variance for the unstable modes as a function of k at four different times: from bottom 
to top t = 30,100,125,150fm/c. The values of o. T, and a are g = 0.3 g eqi T = 4.5MeV, and 
a = 0 . 2 . 


density dependence ( “soft” asymmetry term ). In both cases the density dependence of the 
symmetry energy can be expressed by 

£* 2 ) >-) ((?) ((?2 £?l) 1 


with 


where d 


s{q) = 


2d 


Qeq I T q/ Qeq 

19MeV EnJ, for the “superstiff” case, and 


(24) 


S(q) — di — d 2 g , (25) 

where d\ = 240.9 MeV ■ fm 3 and d 2 = 819.1 MeV ■ fm 6 [22], for the “soft” case. 

The inclusion of the Coulomb interaction presents sizeable effects on the stability con¬ 
ditions of nuclear matter. It gives rise to an overall decrease of the growth rate of density 
fluctuations with a corresponding contraction of the instability region in the (g, T) phase 
diagram [Til. Phi] . Moreover, it can be observed that, when the Coulomb force is included, 
the growth rate vanishes for sufficiently low values of the wave vector k ( k m i n ~ 0.2fm - ) 
0 

In the integrals of Eqs. m and m, which determine the relevant parameters L{t) and 
7 (t) for the distribution P(Z, N,t), we consider only the contributions from the unstable 
modes. To this purpose, we put the weight function f(k) equal to zero for k larger than the 
value beyond which the rate T^ becomes negative. However, to evaluate the total value of 
the covariance matrix, we will consider the sum of the asymptotic value of the contribution 
due to isovector-like fluctuations, Eq. (Cl and the contribution due to the isoscalar-likc 
modes, Eq. 0, that grows exponentially. 

The variance for the unstable fluctuations of the isoscalar density, cr 2 (k) = &fi(k) + 
a 22 (k) + 2of 2 (fc), is displayed in Fig. Q] at four different times. We only report the results 
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obtained with the “superstiff” symmetry term. For the isoscalar fluctuations the “soft” 
asymmetry term gives almost undistinguishable curves. The values chosen for the density 
q = 0.3 Qeq and for the temperature T = 4.5 MeV are in the range expected for the multifrag¬ 
mentation process 34|. For the asymmetry we choose a value of a = 0.2. Figure 0 shows 
that the variance becomes a more and more peaked function about the most unstable mode 
with increasing time. It is worth noticing that the values of the variance of our calculations 
quite well compare with those obtained in Ref. [35J within a different approach including the 
effects of the nucleon-nucleon collisions. This supports the suggestion that the development 
and the growth of the fluctuations are essentially determined by the instabilities of the mean 
held, while the seeds are provided by the thermal agitation of the system. 

We now turn to evaluate fragment isotopic distributions. In order to take into account 
that Z and N are discrete variables we express the probability of finding a correlation domain 
containing Z protons and N neutrons, Y(Z, N,t), through the integral 


rZ pN 

Y(Z,N,t) = / dZ dNP(Z,N,t). 
Jz -1 Jn —1 


(26) 


For large Z and N, Y(Z,N,t) tends to coincide with P(Z, N,t). 

We first consider Eq. (HI to calculate the distribution P(Z, N, t ) and the probability 
Y(Z,N,t). They are determined once the ratio r 0 /L(t) and the parameter y(i) have been 
calculated for given values of g, T and average asymmetry a of the system at the break-up. 
The length L(t) characterizes the decrease of the correlation function with distance. The 
procedure to determine its value has been extensively discussed in Refs. 00 - Here, we 
focus our attention on the calculation of the parameter 7 (t) characterizing the widths of the 
isotopic distributions. 

This can be evaluated by rewriting Eq. CHI) with the assumptions about the weight 
functions introduced in Sec. im 


J dkal 2 (k 1 t)\f(k)\ 2 J dkal 2 (k,t)\f(k)\ 2 
5 dkal 1 (k, t)\f(k)\ 2 f dka\ 2 {k, t) \f(k) | 2 ' 


(27) 


Since the magnitude of the isospin-distillation effect, i.e. the ratio ^i ^k)/af ^k) = 
<j\ 2 {k)/a\ 2 (k), depends on the wave number k, even considering only the contribution of the 
isoscalar-like modes to af J (k) 1 one obtains a non vanishing value of the width 7 . Considering 
also the contribution of isovector-like fluctuations, the width 7 increases, as we will show in 
the following. 

For values of the asymmetry a of nuclear interest, the parameter y(i) turns out to be 
about 10 -3 for both the considered asymmetry terms in the nucleon-nucleon interaction ( 
“soft” and “superstiff” ). As a general trend, the parameter 7 (t) increases with increasing 
asymmetry and density of the decomposing system, and decreases with the time. 

In Fig. 0 we report the isotopic yields of Z — 6 -fragment, calculated according to Eqs. 
m and CHD for two different values of the asymmetry: a = 0.1 and a = 0.2. The used 
values of the parameters T = 4.5 MeV, g = 0.3^ eg and t = 125fm/c, where t is the time 
that the system spends in the instability region, are compatible with the analogous values 
obtained within the SMF approach of Ref. [ 8 j. For the dynamical correlation length we have 
chosen the value L = 1.3 vq. This value corresponds to the effective exponent r e // = 1.65 of 
the power law Y(Z) = YqZ~ T b ff for fragment distribution 0 . In the figure we display the 
results obtained with the “superstiff” asymmetry term and with the “soft” asymmetry term 
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FIG. 2: Calculated isotopic yields of Z = 6 -fragment with the “superstiff” symmetry term (di¬ 
amonds) and the “soft” symmetry term (triangles). The circles represent the results obtained 
neglecting the contribution of isovector-like fluctuations in the “soft” case. The values of g, T, 
L, and t are g = 0.3 g eq , T = 4.5 MeV, L = 1.3 tq, and t = 125fm/c. Top panel: a = 0.1, the 
value of the parameter -y(t) is 7 (f) = 1.02 10~ 3 for the “superstiff” symmetry term, for the “soft” 
symmetry term 7 (t) = 0.69 10 ~ 3 and 7 (t) = 0.3710~ 3 , with and without the contributions from 
the isovector-like fluctuations, respectively. Bottom panel: a = 0.2, the value of the parame¬ 
ter 7 (t) is 7 (t) = 1.62 10~ 3 for the “superstiff” symmetry term, for the “soft” symmetry term 
7 (f) = 0.89 10 ^ 3 and 7 (t) = 0.56 10~ 3 , with and without the contributions from the isovector-like 
fluctuations, respectively. 

of the nucleon-nucleon interaction. Moreover we compare also the relative contribution of 
isoscalar-like and isovector-like fluctuations to the width. 

In the ” superstiff” case isovector-like oscillations are suppressed for the considered values 
of g, T and a, i.e. Eq. (UJ) has only one pole, so the width comes essentially from the disper¬ 
sion of the chemical effect in the isoscalar-like fluctuations (diamonds). In the “soft” case, 
the full calculation is represented by triangles, while the result obtained taking into account 
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only the contribution from the isoscalar-like modes is represented by circles. Comparing 
diamonds and circles, we observe that the “superstiff” asymmetry term gives rise to a wider 
isotopic distribution. This is due to the fact that the “soft” asymmetry term, at the consid¬ 
ered density, is more effective to drive fragments closer to the average asymmetry value, with 
respect to an asymmetry term with a stronger density dependence. Indeed we find that, in 
spite of the competition with Coulomb and surface effects, the isospin distillation mechanism 
does not change much with the wave number k, in the “soft” case. The counterpart in our 
formalism is that in this case the behaviors of the components of the covariance matrix, as 
functions of k, are more similar each other reducing the width of the isotopic distribution. 
However, adding the contributions due to the isovector-like fluctuations, the total width 
obtained in the “soft” case (triangles) becomes closer to the “superstiff” results. It is also 
possible to observe that the contribution of the isovector-like fluctuations to the full width 
is more important at smaller asymmetry. This is because isovector-like fluctuations become 
weaker when increasing the asymmetry of the matter. 

Figure El also shows that the width of the isotopic yields increases with asymmetry. This 
corresponds to the general property that for more neutron-rich systems the density-density 
response function of neutrons is enhanced with respect to that of protons. In addition, 
we can see that the more neutron-rich system (a = 0.2) produces the more neutron-rich 
isotopes, as expected. 

It is worth to remark that both the overall behavior and the widths of the distributions 
of Fig. [21 favourably compare with the corresponding distributions for primary fragments 
calculated within the SMF approach Si]. 


In Fig.[3]we present isotopic distributions obtained using the correlation volume prescrip¬ 
tion (Eq. EEJ), with A = 20. This value corresponds to the average size of intermediate mass 
fragments, as obtained in the considered conditions of density and temperature. As one can 
see by comparing Figs. [21 and |21 results are not very different with the two prescriptions. 

The ratio between isotopic yields observed in two different reactions, Roi(N,Z) = 
Y a2 (N, Z)/Y ai (N, Z), shows a very simple behavior. As a function of Z and N, it can 
well be fitted by an exponential law (the so called isoscaling relationship) [Qjj, |2ll [I2j. In 
addition, the isoscaling relationship has been reproduced by SMF-model calculations also 
for the distributions of primary fragments (23]. This particular feature of the isotopic distri¬ 
butions can represent an effective tool to compare isotopic distributions from systems with 
different N/Z ratios. 

The isotopic ratio R 2 i(N, Z) calculated in our approach, according to Eqs. (fT!l and (l2Hlh 
for two different values of the asymmetry parameter, oq = 0.2 and aq = 0.1, is displayed in 
Figs. 0 and [21 In Fig.Qlwe compare the values of R 2 i(N, Z) as a function of N, obtained with 
the “superstiff” symmetry term and with the “soft” symmetry term. The linear behavior, 
in logaritmic scale, with the same slope for every Z is reproduced in both cases within a 
satisfying approximation. Because of the smaller value of the width parameter y(t), the 
“soft” symmetry term gives a steeper slope with respect to the “superstiff” term. The 
average values of the slope approximative^ are 2.2 ± 0.2 and 1.5 ± 0.15 for the “soft” case 
and the “superstiff” case respectively. 

In Fig. 0 the ratio R 2 i(N, Z) is displayed for two values of the density of the system at 
the break-up. In order to obtain fluctuations of similar magnitude in the two cases, two 
different times the system spends in the instability region are considered. Nevertheless, a 
behavior with a steeper slope is observed in the more unstable case. This is due to a smaller 
value of y(t) in this case, since, for a given charge asymmetry, the response functions of 
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FIG. 3: Isotopic distributions calculated according to the correlation volume prescription (Eq. l2lj) . 
The values of the parameters and the symbols are the same as in Fig El 


protons and of neutrons tend to be more similar with decreasing density. 

We now perform a more quantitative comparison between predictions of our approach 


and results for primary fragments of the SMF-model calculations of Ref. [23]. To this 
purpose we adopt for the average asymmetry of fragments the values predicted by the SMF 
model for semicentral collisions of 112 Sn + 112 Sn and 124 Sn + 124 Sn ja [23[: aq = 0.13 and 
«2 = 0.195, respectively. In both the approaches the same “superstiff” symmetry term 
for the effective interaction is used. Also the values of density g = 0.3 g eq , temperature 
T = 4.5 MeV, and time spent at the break-up t = 125fm/c are chosen according to the 
results of SMF-model calculations. Figure |H1 shows the isotopic ratio R 2 i(N,Z) calculated 
with our approach and the curves obtained in Ref. j23i| by fitting the results of the SMF 
model with an exponential law. We observe a remarkable agreement between the results of 
our nuclear matter calculations and the simulations of the SMF-model. 
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FIG. 4: Isotopic ratio R 2 1 (IV, Z ) = Y a= o^(N, Z)/Y a= o.i(N, Z) calculated with the “superstiff” 
symmetry term (solid lines ) and with the “soft” symmetry term (dashed lines). Lines correspond 
to different values of Z, Z = 3 — 8 from left to right. The values of remaining parameters are the 
same as in Fig 13 



FIG. 5: Same as in Fig.^but using only the “superstiff” symmetry term and for two different values 
of the density: solid lines correspond to g = 0.3 g eq and t = 125fm/c, dashed lines correspond to 
q = 0.4 g eq and t = 150fm/c. 

IV. CONCLUSIONS 

In this article we discuss relevant observables of multifragmentation processes in charge 
asymmetric nuclear matter, such as the isotopic distribution of intermediate-mass frag¬ 
ments, as obtained within the spinodal decomposition scenario, on the basis of an analytical 
approach. Fragmentation happens due to the development of isoscalar-likc unstable modes, 
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FIG. 6: Comparison of calculated isotopic ratio Rm (N. Z) = T a =o.i 95 (-W Z)/l»=o.i 3 (-W Z) 
(diamonds) with the fit for primary fragments of Ref. [23} (solid lines). From left to right 
Z = 3,4, 5,6,7,8. Calculations are done with the “superstiff” symmetry term. The values of 
density g, temperature T, time t, and ratio L/tq are the same as in Fig. [21 


i.e. unstable density oscillations with also a chemical component, leading to the formation 
of more symmetric fragments. We find that the isotopic distributions are peaked at a value 
given by the average distillation effect, while the width is determined by the dispersion of 
the chemical effect among the relevant unstable modes and by isovector-like fluctuations 
present in the matter that undergoes spinodal decomposition. The size of this dispersion is 
mostly due to the competition between symmetry energy effects (that favour the formation 
of symmetric fragments) and the Coulomb repulsion, that acts against the concentration of 
protons in large density domains, expecially for modes with large wavelength. Clearly the 
net result of this competion also depends on the EOS used. Smaller widths are obtained with 
a “soft” symmetry energy term. However, the contribution due to isovector-like fluctuations 
is more important in the “soft” case, indeed in the “superstiff” case isovector oscillations 
are suppressed. Hence finally the isotopic distributions are quite similar when using the two 
parameterizations of the symmetry energy. 

In particular, we find that, when considering two systems with different asymmetry, 
the isotopic (or isotonic) yields obey an approximate isoscaling, with a slope connected to 
the difference betwen the asymmetries of the two systems and to the differences between 
the widths of the isotopic distributions. Hence isoscaling properties can be recovered in 
a dynamical picture. We notice that isoscaling has been found in dynamical simulations 
of heavy ion collisions, such as stochastic mean field and antisymmetrized molecular 


dynamics calculations [3 

The isoscaling parameters are also connected to the properties of the symmetry term 
in the EOS. Indeed we have seen that a stiffer behavior of the symmetry energy term 
yields larger isotopic widths, leading to smaller values of the slope (see Fig. 0 . However, 
as reported in Ref. [|J, we also observe that in collisions of charge asymmetric systems, 
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pre-equilibrium emission is less neutron rich when using a stiffer parametrization of the 
symmetry term (thus leading to more asymmetric fragments), with respect to the “soft” 
case. Therefore, in the isoscaling analysis, there could be a compensation between the 
average asymmetry of fragments (larger in the “stiff” case) and the width of the distribution 
(also larger in the “stiff” case). In fact, for the systems considered in Ref. E3 , similar values 
of the slope are obtained for the two parameterizations considered for the symmetry energy. 

It may also be interesting to notice that the values obtained in our calculations are larger 


than the predictions of statistical multifragmentation models, see Ref. 


22 


Of course this 


picture can be modified by the secondary de-excitation process, that reduces the asymmetry 
of fragments and, consequently, the slopes deduced from isoscaling. Hence the final distribu¬ 
tions can be quite different from the primary ones. A more detailed study, aiming to extract 
information on the primary distributions and on the fragmentation mechanism, would re¬ 
quire the introduction of more sophisticated observables, probably based on an event by 
event analysis, in line with the recent investigations of correlations between intermediate- 
mass fragments 0 - 
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